Pair Phase Fluctuations and the Pseudogap 
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The single-particle density of states and the tunneling conductance are studied for a two- 
dimensional BCS-like Hamiltonian with a d^2_y2-gap and phase fluctuations. The latter are treated 
by a classical Monte Carlo simulation of an XY model. Comparison of our results with recent 
scanning tunneling spectra of Bi-based high- Tc cuprates supports the idea that the pseudogap be- 
havior observed in these experiments can be understood as arising from phase fluctuations of a 
d^2_y2 pairing gap whose amplitude forms on an energy scale set by T*^^ well above the actual 
superconducting transition. 
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Intensive research has focused on the pseudogap 
regime, which is observed in the high-Tc cuprates be- 
low a characteristic temperature that is higher than the 
transition temperature Tc- It occurs in a number of 
different experiments as a suppression of low-frequency 
spectral weight |, |, | |, |, |, |, |. This striking 
pseudogap behavior initiated a variety of proposals as 
to its origin §, |l^, [l6| , since the an- 

swer to this question may be a key ingredient for the 
understanding of high- Tc superconductivity. At present, 
there is no agreement as to which of these proposals is 
correct. In part, this reflects the possibility that there 
may be different pseudogap phenomena operating in dif- 
ferent temperature and doping regimes. In part, this 
is because of the difficulty in determining the experi- 
mental consequences of the various theoretical propos- 
als. In this paper, we focus on the pseudogap phenom- 
ena observed in scanning tunneling spectroscopy mea- 
surements §, § on Bi2Sr2CaCu20s+s (Bi2212) and 
Bi2Sr2CuOQ+s (Bi2201). We provide a detailed numeri- 
cal solution of a minimal model which, however, contains 
the key ideas of the cuprate phase fluctuation scenario: 
that is, we explore the notion that the pseudogap ob- 
served in these experiments arises from phase fluctua- 
tions of the gap 11, 0, H, m, . In this scenario, be- 
low a mean field temperature scale T*^^, a d^2_y2--wave 
gap amplitude is assumed to develop. However, the su- 
perconducting transition is suppressed to a considerably 
lower temperature Tc by phase fluctuations In the 
intermediate temperature regime between T^'^ and 
the phase fluctuations of the gap give rise to pseudogap 
phenomena. 

We will study as a model for phase fluctuations a two- 
dimensional BCS Hamiltonian 



and t denotes an effective nearest-neighbor hopping. The 
{i j) sum is over nearest- neighbor sites of a 2D square 
lattice, and in the second term d connects i to its nearest- 
neighbor sites. In Eq. (Q) one could, of course, add a 
next-near-neighbor hopping t' and a chemical potential 
term. Here, for simplicity and to refrain from further 
approximations, we have set t' and the chemical potential 
equal to zero ||l^. We will assume that below a mean 
field temperature , a dx2-y2-ga.p amplitude forms 
with A ~ 2Tc^^- The detailed temperature dependence 
of A is not central, as we are not interested in the region 
around T^^^ where the pseudogap closes. The important 
point for our calculations is simply that a d2.2_j,2-gap 
amplitude of order 2Tc ^^ in magnitude forms as T drops 
below T*^^ so that 



= A, 



with 



{'Pi 



ipi+s)/2 + TT 



for S in x-direction, 
for S in y-dircction. 
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We then determine the fluctuating phases from a Monte 
Carlo calculation using an effective 2D A"F-free energy 



-£^1 ^cos {(pi 



(4) 



with El adjusted to set the Kosterlitz-Thouless 
transition temperature Tkt equal to some fraction of 
T^ '^- Speciflcally, for the present calculation we will 
set Tkt — T^^^/S. Here, we have the recent scanning 
tunneling results Q for Bi2Sr2CuOe+s in mind, where 
Tc ~ IQK and the pseudogap regime extends to 50 or 
60K, which we take as T*^^. 



2 In principle, the XY action, which determines the fluc- 

H = —t ^ (cJ^Cjcr+c]g.CiCT) — ^ ^(A, ,5(A|^)+A,|^(Ai5)),tuations of the phases, arises from integrating out the 
{ij),cr iS shorter wavelength fermion degrees of freedom including 

(1) those responsible for the local pair amplitude and the in- 
ternal dx2_y2 structure of the pairs. In general this leads 



where cf^. creates an electron of spin a on the z*'* site 
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to a T-dependent quantum action as well as a coupling 
energy Ei , whose temperature dependence is determined 
by the many-body interactions of the microscopic sys- 
tem. There have been various discussions regarding the 
regime over which a classical action is appropriate for the 
cuprates ||l9|, |2^, Here, however, we will proceed phe- 
nomenologically using the classical action, Eq. (^), and 
neglecting the temperature dependence of Ei. Further- 
more, we will use the 2D form of Eq. . One knows that 
for the layered cuprates there is a crossover from 2D to 
3D XY behavior near Tc |22 . Our point of view is that 



away from this crossover regime, a 2D model is certainly 
suitable and on the finite size lattice that we will study, 
the system becomes effectively ordered as T approaches 
Tkt and the correlation length exceeds the lattice size. 
So El will simply be used to set Tkt = Tc- A crucial 
physical point that will be taken into account in our anal- 
ysis is that the basic length scale of the (p-field is larger 
than the Cooper-pair size ^o- Thus, although this is a 
clearly simplified model, we believe that its solution pro- 
vides useful insight into the experimental consequences of 
the phase fluctuation pseudogap scenario. It is the cen- 
tral aim of this paper to verify this by comparison with 
the STM experiments and reproduction of some of their 
characteristic and salient features. 

The calculation of the density of states for an L x L 
periodic lattice now proceeds as follows |2^, |2^. A set 
of phases {ipi} is generated by a Monte Carlo (MC) im- 
portance sampling procedure, in which the probability of 
a given configuration is proportional to exp[—F[ipi\/T) 
with F given by Eq. (^). With {ipi} given, the Hamil- 
tonian of Eq. (|l|) is diagonalized and the single particle 
density of states N{uo, T, {^pi}) is calculated. Further MC 
{tpi} configurations are generated and an average density 
of states N{uj,T) = {N{uj,T, {(fii})) at a given tempera- 
ture is determined. 

As noted above, our point of view is that the XY ac- 
tion, used in the MC simulations, in principle arises from 
integrating out the shorter wavelength fermion degrees of 
freedom up to the scale of the Cooper-pair size, so that 
only the center of mass pair phase fluctuations are im- 
portant. Thus, the scale of the lattice spacing for F[Lpi] 
is set by the pair size coherence length ^ vf/t^^o and 
is of order 3 to 4 times the basic Cu-Cu lattice spacing 
of the fermion Hamiltonian Eq. ([l]). Now the compu- 
tationally intensive part of the calculation is the diago- 
nalization of H and in order to get meaningful results 
as T approaches Tkti we found it necessary to average 
over a large number of Monte Carlo {^Pi\ configurations. 
This requires that some compromise be made with re- 
spect to the lattice size. The results, we will present, 
are for a 32 x 32 Hamiltonian lattice. However, if we 
were to take ~ 4 lattice spacings, this would lead to 
only an 8 X 8 lattice for the ipi simulations. This would 
not allow a sufhcient range for the Kosterlitz-Thouless 
phase coherence length to grow as T approaches Tkt- 



Thus, we have chosen to set A = l.Ot giving ~ 1 so 
that the Lpi simulation can be carried out on the same 
L X L lattice that is used for the diagonalization of H- 
The important physical point is that this procedure effec- 
tively cuts off phase fluctuations on a scale less than the 
Cooper-pair size, ^o- Thus, the phase coherence length 
is always larger than the Cooper-pair size when T is less 
than T*^^. Consequently, our results differ from earlier 
work , which found that the pseudogap regime due to 
fluctuating phases extended only about 20% above Tc, in 
contrast to the Bi tunneling experiments |7j and the 
recent Nernst-effect results . In the work of Ref . , 
parameters were used which set the basic scale of the 
phase correlation length to be much smaller than and, 
therefore, the phase correlation length exceeded only 
in a narrow temperature region set by a fraction of Tkt- 
We believe that this is not the correct phenomenology. 
Results for N{uj,T) are shown in Fig. |l|. For each tem- 
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FIG. 1: Single particle density of states N{ul>) for different 
temperatures T for a 32 x 32 lattice with A = l.Ot and Tkt = 
0-lt- A pseudogap appears below T^'^ ~ 0-5t and coherence 
peaks develop as T approaches Tkt- 

perature we have generated up to 25, 000 independent 
MC {fi} configurations, diagonalized H for each of these 
configurations, and computed {N{Lj,T,{ipi}))- In these 
calculations, as discussed above, we have set A = l.Ot 
corresponding to T*^^ ~ 0.5i and selected Ei so that 
Tkt = 0-lt In order to reduce finite-size effects, 

we employ a very effective scheme recently suggested by 
F. F. Assaad i27|l. 
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For T > , the gap amplitude vanishes and the 
density of states exhibits the usual Van Hove peak at 
u! = Q. For T < T^^^ , the presence of a finite gap ampli- 
tude gives rise to a pseudogap whose size is set by 2A. 
Then, as T approaches Tkt and the XY phase corre- 
lation length rapidly increases, coherence peaks evolve, 
the separation of which is determined by 2A. An impor- 
tant point is that the scale in temperature over which the 
evolution of the coherence peaks occurs, is set by some 
fraction of Tkt which means that it appears suddenly on 
a scale set by T*^^. 

An effective correlation length £,{T), extracted by fit- 
ting an exponential form to the correlation function 

C{e) = (e-*'^'+'e^'^') (5) 

is plotted versus T in Fig. for our 32 x 32 lattice. The 
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FIG. 2: The effective correlation length ^(T) versus T/Tkt 
for the 32 x 32 lattice. Here T^"^ /Tkt ~ 5 so that the 
pseudogap regime which extends from T /Tkt — 1.5 to 5 is 
large compared to the superconducting region which extends 
from to T /Tkt ~ 1. The pronounced increase of ^(T) 
occurs over a narrow temperature region, on a scale set by 
T^'^ , as Tkt is approached. 

rapid onset of £,{T) as Txt is approached is clearly seen. 
It is this sudden increase of ^(T) that is responsible for 
the appearance of the coherence peaks as T approaches 
Tkt- This effect is further enhanced by the 2D to 3D 
crossover that occurs in the actual materials. 

In order to compare these results for N{u},T) with 
scanning tunneling spectra dl/dV, we have calculated 
dI{V, T) / dV using the standard quasi-particle expression 
for the tunneling current. 



dljV^T) 
dV 



NiufJ^d.. (6) 



Here, f{Lo) — (exp(ijj/r) + 1) ^ is the usual Fermi fac- 
tor. Results for dI{V, T)/dV are displayed in Fig. |. The 
effect of the Fermi factors is to provide a thermal smooth- 
ing of the quasi-particle density of states over a region of 
order 2T. This becomes significant at the higher tem- 
peratures and the prominent pseudogap dependence of 




FIG. 3: Tunneling conductance, normalized to its value 
at T^'^ and V = for different temperatures. Solid curves 
are for T = {0.75,1.25,1.75,3.00} Tkt, dashed curves for 
T = {1.00,1.50,2.00} Tkt and T"^ (#|v=o is increasing 
with T). 



N{lu,T) seen in Fig. |i| is smoothed out in dl/dV. In 
Fig. ^, dl/dV results are shown as solid curves for T = 
0.75Tkt (Fig. |a), T = Tkt (Fig. |b) and T - 2Tkt 
(Fig. |c). The dashed curve is for T = T^^^ ~ 5Tkt- 
One sees that the size of the pseudogap scales with the 
spacing between the coherence peaks and evolves contin- 
uously out of the superconducting state. The pseudogap 
persists over a large temperature range measured in units 
of Tkt, becoming smoothed out by the thermal effects 
as T approaches T*^^ and vanishing above T*^^. 

Our numerical results for dI{V,T)/dV are similar to 
recent scanning tunneling measurements of Bi2212 and 
Bi2201 0. Also in these experiments the supercon- 
ducting gap for T < Tkt evolves continuously into the 
pseudogap regime, which extends up to T = T^^ . The 
coherence peaks appear suddenly as Tkt is approached. 
At higher temperatures, the pseudogap fills in rather 
than closing and the temperature range associated with 
the pseudogap regime can be large compared with the 
size of the superconducting regime. 

Summarizing, in order to develop a more quantita- 
tive understanding for the role of phase fluctuations, we 
have provided a numerical solution of a simplified model 
which, nevertheless, contains the key ideas of the cuprate 
phase fluctuation pseudogap scenario. Here the center of 
mass pair-phase fluctuations of a BCS d-wave model were 
determined from a classical 2D XY action by means of 
a Monte Carlo simulation. The resulting tunneling con- 
ductance {dl/dV) reproduces characteristic and salient 
features of recent STM studies of Bi2212 and Bi2201 
suggesting that the pseudogap behavior observed in these 
experiments arises from phase fluctuations of the dx2-y2- 
pairing gap. 
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FIG. 4: Temperature dependence of normalized to its 
value at T^'^ and V — Q. The solid curves are for T = 
0.75Tkt (a), T = Txt (b) and T = 2Tkt (c). The dashed 



curve in all three figures is for T = 



ter (BaCaTeC), the superconducting project KONWHIR 
OOPCV and by the U. S. Department of Energy under 
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